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CHAPTER  1  :  INTRODUCTION 


An  Impact  avalanche  transit  time  (  IMPATT  )  diode  Is  an  attrac¬ 
tive  source  at  millimeter-wave  (  mm-wave  )  frequencies  because  it  Is 
capable  of  producing  a  large  power.  An  IMPATT  diode  Is  a  pn  diode  which 
utilizes  the  avalanche  multiplication  effect  In  the  depletion  layer 
under  the  breakdown  condition  to  generate  mm-wave  energy.  The  size  of 
the  depletion  layer  thickness  of  an  IMPATT  diode  is  determined  by  Its 
operation  frequency.  As  the  frequency  becomes  higher,  a  thinner 
depletion  layer  is  required  for  efficient  IMPATT  operation.  This,  how¬ 
ever,  results  in  a  smaller  negative  resistance  for  the  same  device  area, 
which  makes  impedance  matching  to  the  external  circuit  difficult.  One 
way  to  overcome  this  limitation  is  to  shunt-resonate  the  large  capaci¬ 
tance  of  the  IMPATT  diode.  A  travelling-wave  structure  with  an  appropri¬ 
ate  length  is  one  way  of  realizing  the  shunt  resonance  [1-3].  Fig.  1 
shows  a  very  simple  single-drift  travelling-wave  IMPATT  diode.  The 
structure  is  elongated  in  one  direction  (  the  z  direction  )  and  the  elec¬ 
tromagnetic  wave  travelling  along  this  direction  draws  energy  from  dc 
current  by  the  IMPATT  mechanism  and  Is  therefore  amplified. 

This  type  of  IMPATT  diode  was  first  proposed  by  Midford  et.al., 
[1].  They  actually  constructed  a  silicon  (  Si  )  single-drift  devlced 
and  tested  It  as  an  amplifier.  One  advantage  of  the  travelling-wave 
IMPATT  diode  over  the  lumped  construction  is  that  It  can  be  used  as  an 
oscillating  source  as  well  as  a  two-port  amplifier  depending  on  the 
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Figure  1  Schematic  view  of  the  single-drift 


travelling-wave  IMPATT  diode 
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length  and  the  coupling  condition  to  the  external  circuit.  For  example, 
If  the  device  length  were  one  half  wavelength  at  the  designated  frequency 
and  the  both  ends  of  the  device  were  left  open-circuited,  then  it  would 
oscillate  by  itself  at  the  desired  frequency  without  an  external  reso¬ 
nant  circuit.  This  configuration  makes  the  whole  circuit  simple,  and 
therefore  this  device  is  suitable  for  integrated  circuit.  On  the  other 
hand,  if  the  device  were  connected  to  the  external  circuit  with  good 
impedance  matching,  then  it  would  amplify  the  signal  coming  from  one  end. 
This  amplifier  also  works  in  both  directions,  and  the  structure  is  again 
very  simple. 

The  travelling-wave  IMPATT  diode  has  attractive  features  but  not 
much  work  had  done  on  it  until  the  recent  technological  development  in 
semiconductor  processing.  Molecular  beam  epitaxy  (  MBE  )  is  capable  of 
producing  very  th In  uniformly  doped  semiconductor  layers,  which  is  nec¬ 
essary  for  a  millimeter-wave  IMPATT  diode  [4],  It  is  possible  to  make 
uniform  0.1  ym-thick  layer  over  a  1  mm-long  device  using  MBE  technology, 
and  experimental  work  shows  that  such  a  device  generates  a  useful  amount 
of  power  at  mm-wave  frequencies,  30  -  90  GHz  [5], 

Part  II  comprises  a  theoretical  analysis  of  the  travelling-wave 
IMPATT  diode.  All  possible  losses  caused  by  the  conductors  as  well  as 
the  semiconductors  are  Included  In  the  analysis  since  their  contribution 
at  mm-wave  frequencies  are  not  negligible.  Under  the  small-signal 
assumption,  numerical  results  have  been  obtained  for  SI  single-drift  and 
GaAs  double-drift  travelling-wave  IMPATT  diodes.  Comparison  with  exper¬ 
imental  results  obtained  by  other  authors  shows  good  agreement. 


CHAPTER  2  :  TRAVELLING-WAVE  IMPATT  THEORY 


A  travelling-wave  IMPATT  diode  contains  several  very  thin  semi¬ 
conductor  layers.  This  makes  it  feasible  to  assume  that  the  diode  is 
very  wide  compared  to  the  thickness.  The  structure  to  be  studied,  there¬ 
fore,  reduces  to  a  one-dimensional  multilayered  parallel-plate  waveguide 
containing  an  active  IMPATT  medium.  Fig.  2  shows  the  side  view  of  the 
double-drift  travelling-wave  IMPATT  diode.  Materials  change  only  in  the 
x  direction,  and  the  electromagnetic  wave  is  assumed  to  propagate  in  the 
z  direction.  In  the  y  direction,  both  material  constants  and  field  com¬ 
ponents  are  assumed  to  be  uniform.  Under  an  appropriate  reverse  dc  bias 
condition,  avalanche  breakdown  occurs  mainly  along  the  pn  junction  at 
the  center  of  the  device.  Charge  carriers  travel  in  the  n  and  p  regions 
at  saturation  velocities.  If  a  single-drift  travelling-wave  IMPATT 
diode  Is  considered,  the  two  center  regions  (  n  and  p  )  should  be 
replaced  with  a  single  (  n  )  material.  Since  this  can  be  considered  as  a 
special  case  of  the  double-drift  device,  the  theory  is  developed  for  the 
double-drift  device. 

In  a  multilayered  structure,  the  dominant  propagating  mode  is  a 
TM  mode.  This  mode  has  only  the  y  component  of  the  magnetic  field  and 
the  x  and  z  components  of  the  electric  field.  This  mode  becomes  a  TEM 
mode  If  the  material  contained  between  two  conductors  is  uniform  and  the 
conductors  are  perfect.  However,  the  structure  contains  several  differ¬ 
ent  materials  between  conductors,  and  the  conductors  usually  consist  of 


metal  contact 


Figure 


Side  view  of  the  double-drift 
travel! ing-wave  IMPATT  diode 


a  couple  of  different  metals  In  order  to  form  proper  contacts.  Therefore 
transverse  magnetic  (  TM  )  propagation  occurs  In  the  entire  structure, 
and  establishes  the  required  boundary  conditions. 

The  analytical  process  Is  summarized  as  follows.  A  system  of 
differential  equations  governing  IMPATT  behavior  are  first  solved  to 
determine  the  boundary  of  the  space-charge  region.  For  example,  in  the 
single-drift  diode,  the  space-charge  region  extends  from  inside  of  the 
p+  region  to  the  point  near  the  nn+  boundary  (  Fig.  3  ).  The  wave 
equation  is  then  solved  In  each  layer  Including  this  active  region.  The 
tangential  electric  and  magnetic  fields  must  be  matched  at  every  bounda¬ 
ry.  The  quantity  obtained  from  the  analysis  Is  the  propagation  constant 
of  the  travelling  wave.  This  Is  a  complex  quantity  with  which  the  gain 
and  phase  of  the  wave  In  the  device  can  be  Investigated. 

ANALYSIS  OF  THE  SPACE-CHARGE  REGION 

The  space-charge  region  consists  of  the  avalanche  region  and  the 
drift  region.  In  which  the  dc  electric  field  Is  strong  enough  so  that  the 
carrier  electrons  and  holes  travel  at  saturation  velocities.  In  prac¬ 
tice,  If  the  pn  junction  Is  formed  by  two  uniform  p  and  n  materials,  then 
the  avalanche  and  the  drift  regions  can  not  be  separated.  In  other 
words,  these  two  processes  occur  simultaneously  In  the  space-charge 
region.  In  order  to  physically  separate  these  two  regions.  It  is  neces¬ 
sary  to  control  the  doping  variation  In  the  device  (  for  example,  the 
Read  diode  [6]  ).  In  the  millimeter-wave  IMPATT  diode,  however,  It  Is 


usually  difficult  to  control  the  variation  of  the  doping  level  in  each 
layer  near  the  pn  contact  since  the  entire  structure  is  very  small. 
Therefore,  the  avalanche  and  drift  regions  can  not  physically  be  sepa¬ 
rated.  For  this  reason,  we  have  to  treat  the  space-charge  region  without 
separating  it  Into  two  different  regions. 

The  differential  equations  governing  the  IMPATT  medium  are  as 
follows  [7]: 


2*  * 
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where 
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E  =  electric  field 

-¥ 

J  =  conduction  current 

t 

-► 

J  =  electron  current 

n 

-¥ 

J  =  hole  current 

P 

Ng  =  donor  impurity  concentration 
=  acceptor  impurity  concentration 
n  =  electron  concentration 

p  =  hole  concentration 

a  =  generation  rate  for  electrons 

3  =  generation  rate  for  holes 

q  =  unit  charge 

e  -  permittivity 

p  =  permeability 

The  first  two  equations  are  derived  from  the  Maxwell's  equations  and  the 
remaining  equations  represent  continuity  of  the  current  and  charges. 
The  generation  rate  (  ionization  rate  )  o  and  (5  are  functions  of  the  mag¬ 
nitude  of  the  electric  fields  and  the  relations  are  given  by  [8] 


o 


b 

A  exp  {  -  [ - ]m  } 

|E(x)| 


(  5  ) 


The  values  of  A,  b,  and  m  appearing  In  these  equations  are  given  in 
Table  1  for  Si  and  GaAs.  The  generation  rate  is  the  rate  of  producing 
electron-hole  pairs  which  is  caused  by  avalanche  multiplication. 


To  solve  the  differential  equations  (  1  )  through  (  4  ),  we  have 
to  make  a  few  assumptions.  First  we  separate  each  variable  into  dc  and 
ac  parts.  It  Is  assumed  that  the  ac  signal  is  much  smaller  than  the  dc 
value,  and  we  take  only  the  first-order  sinusoidal  term  in  the  ac  compo¬ 
nent.  Also,  we  assume  a  travelling-wave  form  for  the  ac  term  which  has  a 
propagation  constant  r.  For  example,  the  x  component  of  the  electric 
field  is  written  as 


E 


x 


ejwt-yz 


(  6  ) 


The  dc  bias  voltage  is  applied  only  in  the  vertical  direction  (  x 
direction  ).  Therefore  we  retain  only  the  x  components  of  the  dc  quanti¬ 
ties.  Also,  the  carrier  velocities  in  this  direction  are  limited  to  the 
saturation  velocities  by  the  high  dc  electric  field  in  the  IMPATT  medium. 
The  relation  between  the  x  component  of  the  current  density  and  the  car¬ 
rier  concentration  can  simply  be  described  as: 
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where  vns  and  VpS  are  the  saturation  velocities  of  the  electrons  and  the 
holes,  respectively,  In  the  semiconductor  material . 

Using  these  relations,  the  differential  equations  governing  the 
dc  components  are  written  as  follows: 


dE 


xDC 


dx 


-  (  Nq  -  NA  +  p  -  n  ) 
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dJ 


nxOC 

-  +oJ  n~+3J  nr  =  0 

dx  nxDC  H  pxDC 


dU 


pxDC 

~  °  JnxDC  "  P  JpxDC  =  0 


By  adding  eqs.  (  9  )  and  (  10  ),  we  can  obtain: 


(  8  ) 


(  9  ) 


(  10  ) 


JnxDC  +  JpxDC  =  Jo 


(  11  ) 


where  JQ  Is  a  constant.  Physically,  JQ  represents  the  total  dc  bias  cur¬ 
rent  flowing  In  the  IMPATT  diode.  This  Is  an  arbitrary  constant  and  can 
be  externally  set  to  any  value. 

The  differential  equations  for  the  ac  components  are  obtained  as 


follows: 
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the  following  relations  between  the  z  component  of  the  current  density 
and  the  electric  field: 


J  =  o  E 
nz  n  z 


J  =  o  E 
pz  p  z 


(  16  ) 


where  o  and  o  are  the  conductivities  due  to  the  electrons  and  holes, 
n  p 

respectively.  These  conductivities  can  be  calculated  from  the  distrib¬ 
ution  of  the  dc  current  densities: 


°n  =  *V  =  wn 


o  =  u 
p  *p 


(  17  ) 


where  un  and  yp  are  the  mobilities  of  electrons  and  holes,  respectively. 
Since  the  z  component  of  the  electric  field  can  be  expressed  by  eq. 
(  13  ),  eqs.  (  14  )  and  (  15  )  become: 
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The  ac  equations  to  be  solved  are  eqs.  (  12  ),  (  18  )  and  (  19  ). 


BOUNDARY  CONDITIONS 


In  order  to  solve  the  above  differential  equations,  we  have  to 
specify  all  the  boundary  conditions.  The  boundaries  of  the  space-charge 
region  are  not  yet  known  and  must  be  determined  by  applying  the  dc  bound¬ 
ary  conditions.  The  dc  electric  field  becomes  negligible  outside  of  the 
space-charge  region  since  the  outer  regions  are  much  more  conductive 


than  the  space-charge  region,  which  Is  basically  a  depletion  layer.  At 
the  space-charge  region  boundary  In  the  n  layer,  the  current  consists  of 
only  electron  current.  Also,  at  the  p  layer  boundary,  the  current  con¬ 
sists  of  only  hole  current.  By  applying  these  conditions,  both  bounda¬ 
ries  of  the  space-charge  region  can  be  determined. 

Once  we  have  found  the  boundaries  of  the  space-charge  region,  we 
can  apply  boundary  conditions  to  the  ac  quantities.  The  boundary  condi¬ 
tions  on  the  ac  current  densities  are  the  same  as  in  the  dc  case: 

Jpx  =0  at  the  boundary  in  n  layer 
Jnx  =0  at  the  boundary  in  p  layer 

(  20  ) 

The  tangential  field  components,  Hy  and  Ez>  must  be  matched  to  the  out¬ 
side  regions. 

Outside  of  the  space-charge  region,  we  have  conductive  materials 
which  do  not  contribute  to  producing  rf  energy.  In  such  a  region,  we  can 
write  the  field  components  as 


H  =  A.  e 


+  B,  e 


-j^x 
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where 


rZ  _  „2  .  2 

Ci  -  r  +  « 


and  A.  and  are  constants  determined  from  the  boundary  conditions. 
These  field  components  can  be  matched  at  each  boundary  of  the  metal  and 
semiconductors,  which  leads  to  an  eigenvalue  equation  for  the  propa¬ 
gation  constant.  At  the  space-charge  region  boundaries,  these  field 
quantities  are  transformed  into  field  quantities  used  in  the  ac 
equations  (  12  ),  (  18  )  and  (  19  )  by  using  the  following  relations. 
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obtained  by  Runge-Kutta  integration.  The  second-order  equation  (  12  )  is 
formally  separated  Into  two  first-order  equations,  and  since  the  ac 
equations  are  complex,  they  are  separated  Into  real  and  imaginary  parts 
and  treated  as  different  equations.  Iteration  is  then  performed  until 
all  of  the  above  boundary  conditions  are  satisfied.  The  fourth-  and 
fifth-order  Runge-Kutta  routine  RKF45  In  the  numerical  analysis  library 
at  The  University  of  Texas  was  used  extensively.  This  routine  was  first 
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written  by  H.  A.  Watts  and  L.  F.  Shampine  at  Sandia  Laboratories  and  can 
be  found  In  [9].  Also  the  system  nonlinear  equation  solver,  ZSPOW,  and 
the  complex  root-search  routine,  ZANLYT,  in  IMSL  were  both  used  to  match 
the  boundary  conditions. 


CHAPTER  3  :  SINGLE-DRIFT  TRAVELLING-WAVE  IMPATT  DIODE 


The  theoretical  analysis  was  first  applied  to  the  simple  sin¬ 
gle-drift  travelling-wave  IMPATT  diode  (  n+np+  )  which  had  been  analyzed 
by  Franz  using  a  different  method  [7].  The  material  used  is  Si  and  all 
the  layers  except  the  depletion  region  (  n  layer  )  are  assumed  to  be  per¬ 
fect  conductors.  The  dimensions  and  the  other  parameters  are  shown  in 
Table  2.  In  this  particular  structure,  the  entire  n  region  becomes  the 
space-charge  region  under  operating  conditions.  Therefore  the  boundary 
of  the  space-charge  region  is  already  known.  The  remaining  boundary  con¬ 
ditions  to  be  satisfied  are  the  conditions  on  the  current  densities  and 
the  electric  field:  the  z  component  of  the  electric  field  is  required  to 
be  zero  at  the  n+n  and  np+  boundaries. 

The  calculated  results  are  shown  in  Fig.  4.  In  Fig.  4a,  the  gain 
produced  in  the  diode  is  shown  for  various  dc  current  densities.  The 
upper  half  of  the  figure  indicates  gain  and  the  lower  half  indicates 
loss.  As  the  dc  current  density  increases,  the  gain  also  increases. 
Each  curve  has  a  discontinuity  at  a  certain  frequency.  For  example,  when 
the  dc  current  density  Is  lxlO7  A/m1,  a  large  loss  suddenly  becomes  a 
large  gain  around  6.5  GHz.  The  frequency  of  the  discontinuity  in  the 
gain  corresponds  to  the  avalanche  resonant  frequency  where  the  negative 
resistance  of  the  conventional  lumped  IMPATT  diode  has  a  discontinuity 
[10,11].  The  maximum  gain  constant  is  obtained  just  above  this  disconti¬ 
nuity.  In  practice,  however,  if  the  device  is  operated  at  the  frequency 
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Table  2  Parameters  used  in  the 

single-drift  travelling 

calculation  of  Si 

-wave  IMPATT  diode 

dielectric  constant  (  e„  ) 

r 

11.8 

donor  concentration  (  Ng  ) 

2.2x10* ‘  /m1 

depletion  layer  thickness 

7.0  urn 

saturation  velocity  (  v$  ) 

10s  m/sec 

frequency  GHz 


Figure 


4b  Propagation  constant  versus  frequency  in 
the  single-drift  travelling-wave  IMPATT 
diode  :  Parameter  is  DC  bias  current 


just  above  the  avalanche  frequency,  the  operating  point  moves  back,  and 
forth  between  the  regions  of  gain  and  loss  because  of  the  large-signal 
nature  of  the  actual  device.  Therefore,  the  device  can  not  generate  a 
stable  net  gain.  Hence  the  operating  point  should  be  chosen  well  above 
the  avalanche  resonance  frequency.  This  theoretical  discontinuity 
occurs  because  of  the  small -signal  assumption,  and  therefore  the  theory 
can  not  be  used  to  predict  the  behavior  of  the  device  near  avalanche  res¬ 
onance.  If  the  behavior  of  the  device  near  the  avalanche  resonance  is  to 
be  investigated,  a  large-signal  analysis  needs  to  be  done. 

The  propagation  constant  is  shown  in  Fig.  4b.  The  curves  are 
shown  only  at  frequencies  higher  than  the  avalanche  resonance  frequency. 
The  phase  velocity  of  the  travelling  wave  approaches  the  velocity  of  a 
TEM  wave  in  lossless  Si  material  (  6^=11. 8  )  as  the  frequency  increases. 
Also  the  phase  velocity  Increases  as  avalanche  resonance  is  approached. 

These  results  can  be  compared  to  those  obtained  theoretically  by 
Franz  [7],  The  theoretical  values  of  gain  are  much  smaller  than  Franz's 
while  the  behavior  of  the  propagation  constants  agrees  qualitatively. 
The  difference  is  due  to  different  simplifying  assumptions.  Fig.  5 
shows  the  dc  solution  of  the  present  device.  Fig.  5b  indicates  that  ava¬ 
lanche  multiplication  occurs  in  more  than  half  of  the  entire  depletion 
region.  This  means  that  the  assumption  made  by  Franz  that  the  avalanche 
region  is  infinitesimally  thin  and  the  entire  space-charge  region  can  be 
treated  as  a  drift  region  should  not  be  used.  Actually  carrier  drift  and 
avalanche  multiplication  occur  simultaneously.  It  is  not  necessary  to 
distinguish  these  two  regions  in  this  analysis.  The  effect  of  the  loss 
caused  by  the  inactive  regions  will  be  discussed  in  the  next  chapter. 


distance  pm 


Figure  5a  Solution  of  the  DC  equation  :  DC 

electric  field  is  plotted  as  a  function 
of  the  distance  from  the  n+n  junction. 
Parameter  is  DC  bias  current  density. 


distance  pm 


Figure  5b  Solution  of  the  DC  equation  :  DC  hole 

current  density  is  plotted  as  a  function 
of  the  distance  from  the  n+n  junction. 
Parameter  is  DC  bias  current  density. 


CHAPTER  4  :  DOUBLE-DRIFT  TRAVELLING-WAVE  IMPATT  DIODE 


A  double-drift  IMPATT  diode  Is  potentially  capable  of  producing 
more  power  than  a  single-drift  diode  [12-14].  The  depletion  region  of 
the  double-drift  diode  consists  of  two  layers  (  p  and  n  ).  The  dc  elec¬ 
tric  field  Is  maximum  at  the  center  junction  (  pn  junction  )  and  there¬ 
fore  avalanche  multiplication  primarily  occurs  in  this  region. 

The  parameters  chosen  for  the  numerical  calculations  are  shown 
in  Table  3.  The  structure  consists  of  8  layers  and  Is  very  similar  to 
the  one  used  In  the  experiment  of  Bayraktaroglu  et.al.  [5].  The  Inactive 
layers  are  treated  as  lossy  dielectric  layers  and  their  effect  is  taken 
into  account. 

The  computed  results  for  gain  and  propagation  constants  are  pre¬ 
sented  in  Fig.  6.  As  frequency  decreases  and  approaches  avalanche  reso¬ 
nance,  the  gain  increases  and  the  phase  velocity  of  the  wave  increases. 
In  these  figures,  dotted  lines  Indicate  the  results  for  the  special  case 
that  the  depletion  region  Is  replaced  by  a  lossless  material,  where 
slow-wave  propagation  becomes  evident.  Each  curve  asymptotically 
approaches  the  dotted  line  at  high  frequencies,  since  the  device  becomes 
Inefficient  there  and  deviates  from  IMPATT  operation. 

Fig.  7  compares  the  present  theory  and  the  experimental  results 
obtained  by  [5].  The  device  was  used  as  an  oscillator  with  one  end  open 
and  the  other  terminated  In  a  short.  In  the  figure,  white  circles  (o) 
show  the  actual  device  lengths  and  the  oscillation  frequency  obtained  in 


Table  3  Parameters  used  In  the  calculation  of  GaAs 

double-drift  travelling-wave  IMPATT  diode 


material  thickness  (pm)  doping  level  (/m*)  conductivity  (S/m) 


Au 

7.0 

- 

4.3x10* 

Ti 

0.1 

- 

u> 

O 

rH 

X 

CO 

r—i 

n+  GaAs 


.2 


5.0x10** 


6.8x10s 


frequency 


Propagation  constant  versus  frequency  in  the 
double-drift  travelling-wave  IMPATT  diode  : 
Parameter  is  DC  bias  current  density. 

-  :  Active  region  is  replaced  by 

lossless  GaAs  material. 


oscillation  free 


Figure  7  Comparison  between  experiment  and  theory  :  Solid 
lines  represent  the  region  where  net  gain  is 
produced  in  the  device. 


—  :  Active  region  is  replaced  by  loss¬ 
less  GaAs  material. 
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the  experiment,  and  black  dots  (•)  show  one-third  of  the  actual  device 
length  (  i.e.,  the  device  Is  considered  to  be  oscillating  at 
three-quarter  wavelength  ).  The  three  solid  curves  indicate  the  theore¬ 
tical  quarter  wavelength  where  net  gain  is  produced  in  the  device.  They 
are  In  good  agreement  except  for  two  black  dots  where  higher  order  reso¬ 
nance  (  3X/4  )  takes  place.  In  these  cases,  the  actual  terminating  con¬ 
ditions  might  have  a  more  dominant  effect  on  the  oscillation  frequencies 
(  i.e.,  not  quite  open  circuit  ). 

Fig.  8  shows  how  gain  changes  with  respect  to  the  dc  current.  As 
the  dc  current  increases,  the  gain  also  Increases  but  suddenly  goes  to  a 
large  loss  below  a  critical  bias  point.  This  fact  qualitatively  agrees 
with  the  results  of  power  measurements  for  conventional  IMPATT  diodes 
[♦]• 

The  oscillation  frequency  of  the  double-drift  travelling-wave 
IMPATT  diode  Is  plotted  with  respect  to  the  dc  bias  current  in  Fig.  9. 
The  device  length  is  fixed  at  30  mils.  Circles  are  the  experimental 
results  obtained  by  Bayraktaroglu  [15].  Our  theoretical  curve  closely 
follows  them.  In  their  experiment,  a  mode  jump  was  observed  around 
Jq=8x107  A/m*,  although  the  theory  shows  continuous  transition.  In 
order  to  investigate  this  phenomenon,  it  Is  really  necessary  to  include 
the  effect  of  the  coupling  to  the  external  circuits.  In  the  figure,  the 
dotted  curve  shows  the  oscillation  frequency  If  the  propagation  constant 
Is  equal  to  that  of  TEM  wave  In  a  lossless  GaAs  material  (  er=13.0  ). 
Therefore,  the  theoretical  and  experimental  results  show  the  fast  wave 
nature  of  the  travelling-wave  IMPATT  diode. 


CHAPTER  5  :  PERIODIC  IMPATT  OSCILLATOR 


So  far,  the  discussion  pertains  only  to  the  uniformly  distrib¬ 
uted  IMPATT  diode.  It  was  shown  that  such  a  diode  produces  gain  in  the 
electromagnetic  wave  travelling  along  the  pn  junction.  Such  a  diode  is 
required  to  have  extremely  uniform  semiconductor  layers  in  order  to  have 
the  same  breakdown  voltage  over  the  diode  area. 

A  travelling-wave  IMPATT  oscillator  can  also  be  realized  by 
using  a  conventional  TEM  parallel-plate  waveguide  loaded  periodically 
with  lumped  IMPATT  diodes.  This  structure  is  shown  in  Fig.  10.  Since 
each  diode  has  a  small  area,  it  is  easy  to  fabricate.  Also,  the  struc¬ 
ture  Includes  the  resonant  circuit  if  appropriately  cut  and  terminated. 
Therefore,  we  again  do  not  need  any  external  resonant  circuit,  and  the 
entire  system  is  kept  simple. 

In  this  chapter,  a  simple  analysis  of  the  periodic  IMPATT  oscil¬ 
lator  is  presented.  The  structure  Is  basically  modeled  by  a  transmission 
line  loaded  with  IMPATT  Impedances.  The  effect  of  the  supporting  metal 
posts  are  also  Included  by  calculating  the  inductance  of  the  post.  Some 
preliminary  results  are  shown  in  the  last  section  of  this  chapter. 


CALCULATION  OF  THE  POST  IMPEDANCE 


The  parallel-plate  waveguide  loaded  periodically  with  IMPATT 
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diodes  Is  modeled  by  an  equivalent  transmission  line  loaded  with  Impe¬ 


dance  Zp.  The  Impedance  Zp  Is  a  sum  of  the  ac  IMPATT  Impedance  and  the 
metal  post  reactance.  In  this  work,  we  treat  the  IMPATT  diode  as  a 
lumped  element,  and  the  electric  field  component  and  current  density  are 
assumed  to  be  oriented  In  the  x-dl recti  on.  The  governing  dc  equations 
are  the  same  as  those  developed  In  Chapter  2,  l.e.,  (  8  ),  (  9)  and 
(  10  ).  The  ac  equations  have  simpler  forms  as  follows: 


dE  1  1  1  ~  J  jwe  ~ 

_  =  -  {  (  —  +  —  )  J - +  —  E  } 

dx  '  V  vns  vns  vns 

(  23  ) 

jw  ~ 

(  3  -  a  -  -  )  Jp  +  aJac 

V 

+  (  ^nOC®'  +  JpDC*'  "  juca  )  E 

(  24  ) 

These  equations  are  solved  In  a  similar  way  to  the  previous  analysis  of 
the  travelling-wave  IMPATT  diode.  The  ac  electric  field  Is  then  inte¬ 
grated  over  the  entire  IMPATT  region  to  give  the  ac  voltage,  from  which 
we  can  obtain  the  ac  Impedance  of  the  diode. 

The  metal-post  reactance  needs  to  be  calculated.  Since  the 
thickness  of  the  IMPATT  diode  Is  much  smaller  than  the  length  of  the  sup- 


% 

dx 


porting  metal  post,  we  can  not  neglect  the  effect  of  the  reactance  of  the 
post.  The  metal  post  is  assumed  to  be  uniform  in  diameter.  We  use  the 
conventional  mode-matching  technique  to  calculate  the  first-order  induc¬ 
tance  of  the  post  [16].  The  cross-sectional  view  of  the  waveguide  model 
Is  shown  In  Fig.  11.  The  post  Is  assumed  to  be  located  at  2=0.  A  magnet¬ 
ic  wall  is  placed  on  both  sides  and  the  field  around  the  metal  post  can 
be  expanded  as  follows: 

~  rrmx  -2f  1  z  | 

E  =  I  A  cos —  e  m 
y  m=0  m  a 

(  25  ) 


1  «•  mirx  -If  |z| 

H  =  -  l  -  r  Am  cos —  e  sgn(z) 

x  j«y  m=0  m  m  a 

(  26  ) 


where 


The  assumption  of  continuity  of  H  at  z=0  gives: 


circular  post 


Figure  11  Cross-sectional  view  of  the  analytical  model 

for  the  calculation  of  the  metal -post  reactance 


—  £  Am  jr  cos — 
jwy  m=0  a 


-  =  I  6(x-D) 


(  27  ) 


where  I  Is  the  current  flowing  on  the  metal  post  and  6(x)  is  the  Dirac 

delta  function.  Using  Fourier  analysis,  we  can  obtain  A  and  hence: 

m 


-  jwp  mirD  mirx  —  JT  |  z  | 

E  =  l - I  cos —  cos —  e 

y  m=0  T  a  a  a 

m 


(  28  ) 


We  now  assume  that  the  Incident  wave  is: 


E1  =  « 


The  tangential  electric  field  on  the  surface  of  the  metal  post  must  equal 
zero.  We  only  choose  one  point  (  x=D+r,  z=0  )  to  do  this: 


Hence  we  can  obtain  the  expression  of  the  current  I 
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«  r  «•  mirr  2mnD  1 

=  —  £  {  cos —  +  cos -  }  /  rl 

jwy  L  m=0  a  a  m  -* 


Using  this  result,  the  impedance  of  the  post  jXQ  is 


(  30  ) 


k  a  irr  nD 

X  =  -  {  —  log  (  4  sin —  sin —  ) 


4  ir 


2a  a 


-  2mirD  1  a 

+  £  (  1  +  cos -  )  ( - —  )  > 


a  y  mu 
m 


If  the  post  is  located  at  the  center  (  D=a/2  ),  then 


k  a  2nr  «•  la 

XQ  =  -  {  -  -  log( — )  +  £  2  ( - )  > 

4  ir  a  m=2  y_  mir 


(  31  ) 


(  32  ) 


These  Impedances  are  normalized  values  with  respect  to  the  character¬ 
istic  impedance  of  the  parallel-plate  waveguide.  The  total  post  impe¬ 
dance,  Zp,  is  the  sum  of  the  ac  IMPATT  impedance  and  this  post  reactance. 


RESULTS 


Once  the  value  of  the  loaded  Impedance  Zp  Is  known,  the  charac¬ 
teristics  of  this  structure  can  be  found  using  the  theory  of  periodic 
structures.  However,  a  different  approach  is  used  here.  In  practice,  it 
is  not  necessary  to  have  many  diodes  and  it  is  also  desirable  to  keep  the 
total  line  length  short  to  ensure  operation  at  one  frequency.  That  is, 
the  separation  of  the  resonant  frequencies  becomes  large  as  the  length 
decreases.  As  Fig.  12  shows,  simple  admittance  transforms  are  used  to 
analyze  the  resulting  structure.  For  example, 


yL  *  tanh  Ydout 
1  *  Yl  tanh  Ydout 


Y  +  tanh  Yd. 

—  n  in 

Yin  =  — - 

1  +  Y„  tanh  Yd. 
n  in 

(  33  ) 

The  device  should  oscillate  at  the  frequency  where  the  imaginary  part  of 

the  input  admittance  Y.  becomes  zero. 

in 


As  an  example,  a  parallel-plate  waveguide  loaded  with  four 

IMPATT  diodes  is  chosen.  Its  dimensions  and  the  IMPATT  parameters  are 

shown  in  Table  4.  The  waveguide  is  filled  with  polyimide  which  has  a 

dielectric  constant  of  e  =3.4,  and  the  device  ends  are  assumed  to  be 

r 

open-circuited.  The  calculated  input  admittance  is  shown  in  Fig.  13. 
Oscillation  occurs  at  50  GHz  in  the  figure  where  the  imaginary  part  of 
Yin  becomes  zero  and  the  real  part  negative  maximum. 
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Table  4  Parameters  of  periodic  IMPATT  oscillator 


length  of  parallel-plate  waveguide 
width  of  parallel-plate  waveguide 
distance  between  diodes 


IMPATT  parameters 
p+  GaAs  layer  (  Np=5.0xl02 */m3  ) 
p  GaAs  layer  (  Nq=1.5x102 Vm*  ) 
n  GaAs  layer  (  N^=1.5xl023/m3  ) 
n+  GaAs  layer  (  N^=5.0xl02Vm3  ) 


50  mils 
10  mils 


10  mils 


0.1  pm 
0.35  pm 
0.35  pm 

0.1  pm 


.-\0 
.  vi 
'■  1 
i  _  ... 


Figure  13  Computed  input  admittance  of  the  periodic 
IMPATT  oscillator  :  Resonant  frequency  is 
about  50  GHz. 


CHAPTER  6  :  CONCLUSIONS 


Part  II  covered  the  theoretical  discussion  of  the  travel- 
ling-wave  IMPATT  diode.  Small-signal  theory  governing  the  behavior  of 
an  electromagnetic  wave  travelling  in  such  a  structure  was  developed. 
The  theory  showed  the  conditions  where  the  device  produces  gain.  Using 
the  theory,  one  can  predict  the  performance  of  the  travelling-wave 
IMPATT  diode.  Experimental  results  obtained  elsewhere  confirm  the  theo¬ 
ry,  although  the  device  behavior  near  avalanche  resonance  could  not  be 
predicted  well . 

The  periodic  IMPATT  oscillator  was  also  discussed  as  another 
form  of  the  travelling-wave  structure.  A  simple  theory  predicted  the 
oscillation  frequency. 
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APPENDIX  A  :  PROGRAM  LISTING 


Appendix  A  contains  the  FORTRAN  program  code  for  the  analysis  of 
the  propagation  characteristics  of  the  travelling-wave  double-drift 
IMPATT  diode.  The  subprogram  RKF45  is  found  in  the  NATEXT  library,  and 
ZSPOW  and  ZANLYT  are  found  in  the  IMSL  library  in  The  University  of  Texas 
Numerical  Analysis  Library. 


uuuuuouuuouuuouuuuuuuuuuuuuuuuuuuuuuuuuuuu  uuu 
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PROCRAM  IMPATT  <TTY. OUTPUT, TAPE!. TAPE5-TTY. TAPEA-TTY) 


COMPUTATION  OF  THE  COMPLEX  PROPAGATION  CONSTANT 
OF  THE  TRAVELL INC -WAVE  DOUBLE-DRIFT  IMPATT  DIODE 


STRUCTURE  : 

THE  CROSS  SECTION  OF  THE  STRUCTURE  IS  SHOWN 
BELOW: 


GOLD 

T1 

TITANIUM 

T2 

N+  GAAS 

TN 

N  GAAS 

TP 

u> 

< 

s 

CL 

T3 

P+  GAAS 

T4 

TITANIUM 

GOLD 

T1  -  T4  .  , .  THICKNESS  OF  EACH  LAYER 

(VALUES  ARE  ASSIGNED  IN  THE  PROCRAM.  ) 
IMPURITY  DENSITY  OF  EACH  LAVER  IS  GIVEN  BY 
THE  FUNCTION  OENS(X). 

GOLD  LAYERS  ARE  ASSUMED  TO  BE  SEMI -INFINITELY 
THICK. 


LIBRARIES  :  NALIB4  <RKF45> 

IMSLIBF  (ZSPOU, ZANLYT> 


COMMON 

/BASE /PI,  DCJ.  G,  EPSO,  AMUO.  OMEGA 

COMMON 

/BA3EM/VP, VN, CMUM, GMUP.  EPSG 

COMMON 

/ACTIVE/TN,  TP 

COMMON 

/BOUND/ XI, XF 

COMMON 

/LAVER/ZEPSO. ZEPS, T 

COMMON 

/PROP/ZGAM, ZD 

COMMON 

/ACEXI/ZEXI 

COMPLEX 

ZEPSO,  ZEPS 

COMPLEX 

ZJ.  ZGAM.  ZEXI,  ZD.  ZDET,  CX 

DIMENSION 

ZEPS ( 6 ) , CX (2) .  INFER ( 2 ) 

DIMENSION 

X (2) , F( 2 ) , PAR (2) , WORK (21 > ,  S(6) ,  T(6> 

EXTERNAL 

FCNJ ,  ZDET 

BASIC  CONSTANTS. 

ZJ  «  10.0,1.0) 

PI  -  3. 1415924S059 
DBN  -  8. 6859E-3 


UUU  UUU  UUUU  U  u  u  uuuu 
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QMUN  «  0.  83 
GMUP  -  0  04 
0  -  1. 602E-19 
EFSO  -  e.  B54E-12 
EPSC  -  13.  0 
AMUO  ■  4.  0*PI/1.E7 
VP  -  6.  0E4 
VN  »  6.  0E4 

THICKNESS  OF  EACH  LAVER. 

T1  *  0.  IE-6 
T2  -  0.  2E-6 
T3  ■  0.  2E-6 
T4  -  0.  IE-6 
TN  »  0.  3E-6 
TP  -  0.  3E-6 

CONDUCTIVITIES. 

SCOLD  -  4.  3E7 
SCI)  **  1.  8E6 

SC2)  -  Q*GMUN*DENS( -TN/2.  0) 

SC3)  -  Q*GMUN*DENS(TN/2.  0) 

S<4>  -  -G*CMUP*DENSCTN+TP/2.  0) 

S(S)  -  -G*GMUP*D£NS ( TN+TP+TP/2.  0) 

S(6)  -  1.  8E6 

START  DC  CALCULATION  WHICH  DETERMINES 
THE  BOUNDARIES  OF  THE  ACTIVE  REGION. 

1000  WRITE  <6.  100) 

100  FORMAT  (  /3X. 'DC  CURRENT  <A/M2> '  > 

READ  (9.  •  )  DC J 

NSIG  -  8 
N  -  2 
ITMAX  »  90 
X(l)  -  0.  03*TN 
X(2>  -  0.  95*<TN+TP> 

CALL  ZSPOW  (FCNl,  NSIG.  N.  ITMAX.  PAR.  X.  FNORM. 
1  WORK.  IER) 


XI  -  X<1) 

XF  -  X (2) 

WRITE  (6.200)  XI,  XF,  IER 

200  FORMAT  <  /3X.  'XI- '.  E12  6,  3X,  'XF- E12.  6.  3X,  '  IER-  M3  ) 

DETERMINE  BOUNDARIES  FOR  THE  CALCULATION 
OF  THE  PROPAGATION  CONSTANT. 

T(  1 )  -  T1 
T<2>  -  T2 
T<3)  -  XI 

IF  C  XI  .  QE.  0.0)  GO  TO  10 
T<2)  -  T2+XI 
T<3)  -  0.  0 
10  T(4i  ■  TN+TP-XF 
T<3)  ■  T3 

IF  l  XF  LT  TN*TP  l  CO  TO  20 
T<4)  ■  0.  0 


T(5>  -  T3+TN+TP-XF 
20  T(6)  -  T4 

START  AC  CALCULATION  WHICH  DETERMINES 
THE  COMPLEX  PROPAGATION  CONSTANT. 

2000  WRITE  (6.110) 

110  FORMAT  <  /5X.  'FREQUENCY  <CHZ>'  > 

READ  (5.  *  )  FREQ 

OMEGA  -  2. 0£9*PI*FREG 
AKO  »  OMEGA/3.  OES 
ZEPSO  -  EPSO-Z J*SGGLD/OMEGA 
ZEPS( 1 )  -  EPSO-Z J*S<1) /OMEGA 
ZEPS(2>  »  EPS0*EPSG-ZJ*S<2) /OMEGA 
ZEPSO)  -  EPS0*£PSG-ZJ*S(3> /OMEGA 
ZEPS<4)  -  EPSO*£FSG-Z J*S<4 ) /OMEGA 
ZEPS(S)  «  EPS0*EPSG-ZJ*S(3> /OMEGA 
ZEPS16)  -  EPSO-Z J*S(6) /OMEGA 

NSIC  *  3 
ITMAX  -  25 
30  WRITE  (6.230) 

230  FORMAT  (  /SX, 'INITIAL  GUESS  :  ALPHA.  BETA.  1/2 

1  /21X.  '1  .  .  .  NEPER/M.  /M' 

2  /21X,  '2  ..  DB/MM,  BETA/KO'  > 

READ  <5.  •  )  ALPHA1, SETA1.  KK 
IF  (  KK  .  EQ.  1  )  CO  TO  40 

IF  (  KK  NE  2  )  GO  TO  30 
ALPHA 1  -  ALPHA1/DBN 
BETA1  -  BETAKAKO 
40  CX<1)  -  ALPHA 1  +  ZJ*BETA1 
ZEXI  =  (1.0,  1.  0) 

CALL  ZANLYT  (2DET .  0  0.  NSIG.  0,  1,  1.  CX.  ITMAX. 

1  INFER. IER) 

IF  (  AIMAG  (CX  ( 1 )  )  .LT.  0.0  >  CX(1>  -  -CX(1> 

ALPHA 1  »  REAL (CX ( 1 ) ) 

BETA1  =  AIMAG(CX  < 1 ) ) 

ALPHA2  »  ALPHA1*DBN 
BETA2  -  BETA 1 /AKO 

WRITE  (6.240)  DCJ.  XI.  XF.  FREQ, 

1  ALPHA1. ALPHA2. BETA1.  BETAZ.  IER 

240  FORMAT  (  //SX.  'DC  CURRENT  -  '.EIO.  4,  '  A/M' 

1  /8X,  '(  XI  -  '•  E10.  4,  '.  XF  -  EIO.  4.  '  )  ' 

2  /SX, 'FREQUENCY  -  ',F6. 2. '  GHZ' 

3  //5X.  'ATTEN  CONST.  -  '.EIO.  4.'  NEPER/M,  '. 

4  EIO.  4,  '  DB/MM '/SX.  'PROP  CONST.  *  '.EIO.  4, 

5  '  /M,  BETA/KO  - ' .  F10  6.  SX.  '  IER  -'.13  ) 

WRITE  (1.240)  DCJ. XI. XF, FREG, 

1  ALPHA 1,  ALPHAS.  BETA 1 ,  BETA2,  I ER 

SO  WRITE  (6.230) 

2S0  FORMAT  (  /SX,  '1  ...  CHANGE  DC  CURRENT' 

1  /SX, '2  ..  CHANGE  FREQUENCY' 

2  /SX, '3  ...  CALCULATE  FIELD  COMPONENTS' 

3  /SX, '4  ...  STOP'  ) 

READ  <S.  *  )  ICTRL 

IF  <  ICTRL  . EQ  l  )  GO  TO  1000 


oouuu  ooo  g  u  u  u  u  uuuuu  uuu 
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IF  (  I CTRL  .  EQ.  2  )  SO  TO  2000 
IF  <  ICTRL  .  NE.  3  >  STOP 
C 

CALL  FLDPRT  (50) 

C 

CO  TO  50 
C 

END 


SUBROUTINE  FCN1  (X.F, N. PAR) 


SOLVE  DC  EQUATIONS 


DIMENSION 

COMMON 

COMMON 

COMMON 

EXTERNAL 


X<2).  F<2).  PAR<2>.  Y<2>.  W0RK<63>.  IMORK(S) 

/RK JUNK /WORK. IWORK 

/BASE/PI,  DC J.  Q.  EPSO.  AMUO,  OMEGA 

/BA3EM/VP,  VN.  CMUN,  CMUP,  EPSC 

DCEQN 


NEON  -  2 
I FLAG  -  1 
RELERR  -  I. OE-8 
ABSERR  «  0.  0 


XI  «  XU) 
X2  -  X (2) 
Y<1 )  -  0.  O 
Y<2>  -  DCJ 


CALL  RKF43  (DCEQN. NEON,  Y. XI, X2.  RELERR, ABSERR. 

1  IFLAC.  WORK.  IWORK) 

IF  (  IFLAC  NE  2  )  CO  TO  1000 
F<1)  -  YU) 

F(2>  -  Y(2) 

WRITE  (6.200)  XU).  X  (2) .  FU  >,  F(2) 

200  FORMAT  <  3X,  'XI=',E9.  3.  '  XF«',E9.  3,  '  DCEX-'. 

1  E9.  3,  '  DCJN-',E9.  3  ) 

RETURN 

1000  WRITE  (6. 1010)  IFLAC 

1010  FORMAT  <  /5X, 'ERROR  (FCN1-RKF45)  ...  IFLAC-', 13  ) 
RETURN 
END 


SUBROUTINE  DCEON(  X,  Y.  YPRIME) 

DC  EQUATIONS 

DIMENSION  Y(2K  YPRIME<2) 

COMMON  /BASE /PI, DCJ,  Q,  EPSO.  AMUO, OMEGA 

COMMON  /BASEM/VP. VN. CMUN,  CMUP.  EPSC 


nnn  o  o  n  n  o 
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DCJN  -  Y<2> 
C 


YPRIME(l)  -  ( Q*rENS< X  )  *DC J/VP-DC JN*< 1  /VP+1  /VN>> 
t  /EPSO/EPSG 

YPRIciE;2)  -  -CENA (DCEX  >*DC JN-QENB (DCEX  >*(DCJ-DC JN) 

C 

RETURN 

END 

C 

C 

C 

C 

C 

COMPLEX  FUNCTION  ZDET  ( Z6AMMA ) 

C 

C  SOLVE  FOR  THE  PROPAGATION  CONSTANT 

C 

COMPLEX  ZCAMMA,  ZGAM,  ZEXI,  ZD,  ZJ 
COMMON  /PROF /ZGAM. ZD 

COMMON  /ACEXI/ZEXI 

EXTERNAL  FCN2 

C 

ZGAM  -  ZGAMMA 
EPS  *  0.  0 
NSIG  ■  5 
ITMAX  *  29 
C 


CALL  ZCOMP  <FCN2,  ZEXr.EPS. NSIG.  ITMAX.  ICOND) 
ZDET  *  ZD 

WRITE  (A. 200)  ZGAM. ZDET. ICOND 
200  FORMAT  <  9X.  'ZCAM  ■ ' .  E9.  3.  X.  E9.  3.  '  ZDET 
1  E9.  3.  X,  E9.  3.  '  I  COND=»  ',11  > 

RETURN 

END 


SUBROUTINE  FCN2  < ZJN, ZEXI ) 


C 


C 

c 

c 


MATCH  BOUNDARY  CONDITIONS. 


DIMENSION 

DIMENSION 

COMPLEX 

COMPLEX 

COMMON 

COMMON 

COMMON 

CCMMON 

COMMON 

COMMON 

EXTERNAL 


Y<  10) .  WORK  (63).  I  WORM  3) 

T<SS. ZEPS<6) 

ZEPSO.  ZEPS,  ZGAM,  ZEXI,  ZD,  ZKX,  ZJ 

ZJP,  ZJN,  ZEX,  ZEXP.  ZHY,  ZEZ,  ZZl,  ZZ2.  ZA.  ZB 

/BASE/PI,  DCJ,  0.  EPSO.  AMUO,  OMEGA 

/BA5EH/VP, VN,  GMUN,  OMUP,  EPS6 

/RK  JUNK /WORK,  IWORK 

/BOUND/XI, XF 

/PROP /ZGAM. ZD 

/LAYER /ZEPSO,  ZEPS,  T 

ACEGN 


ZJ  -  <0  0. 1  0) 

ZEX  -  ZEXI 

UPPER  GOLD  REGION 


ZKX  m  CSGRT ( Z GAM# ZGAM+OMEGA(*OMEGA* ZEPSO* AMUO) 


I 


I 


non  non 


IF  <  AI MACi ZKX  >  .  GT.  0.0  >  ZKX  -  -ZKX 
ZHY  -  (1.0. 0  0) 

ZEZ  -  ZKX /OMEGA/ ZEPSO 

TITANIUM.  N+  AND  N  REGIONS 

DO  30  1*1. 3 

IF  '  T<  I  >  EG.  0  0  )  CO  TO  30 

ZKX  -  CSGRT ( ZGaM*ZGAM+OMEGA*OMECA*ZEPS ( I >*AMUO) 

ZA  ■  0.  3*( ZHY+OMEGA*ZEPS( I >* ZEZ /ZKX  > 

ZB  *  0.  5*(ZHY-0MECA*ZEPS(I)*ZEZ/ZKX> 

ZHV  *  ZA*CEXP(ZJ»ZKX*T(I>>  +  ZB*CEXP(-ZJ*ZKX*T( I > > 
ZEZ  -  (ZA*CEXP(ZJ*ZKX*T(1))  -  ZB*CEXP(-ZJ*7KX*T( I > ) > 
1  *ZKX/OMEGA/ZEPS(I> 

30  CONTINUE 

ACTIVE  REGION 

NEGN  -  10 
I FLAG  »  1 
RELERR  *  1. 0E-8 
ABSERR  »  1.  OE-12 
XI  -  XI 
X2  =*  XF 

ZJN  *  ZGAM*ZHY-ZJ*OM£GA*EPSOi*EPSG*ZEX 

ZEXP  -  ZCAM*Z£Z-Z  JN/VN/EPSO-'EPSG 

Y(  1  >  *  0.  0 

Y(2)  *  DCJ 

Y<3>  -  0.  O 

Y<4)  ■  0.  O 

Y< 5 )  ■  REAL (ZJN) 

Y( A )  -  AIMAC(ZJN) 

Y( 7 )  -  REAL(ZEX) 

Y(8)  ■  AIMAC(ZEX) 

Y(9 )  -  REAL ( ZEXP ) 

Y ( 10)  -  AIMAG(ZEXP) 

CALL  RKF45  ( ACEQN. NEGN. Y.  XI, X2,  RELERR.  ABSERR. 

1  I FLAG. WORK, I WORK) 

IF  (  I FLAG  .  NE  .  2  >  GO  TO  1000 
ZJP  ■  Y(3)  ♦  ZJ*Y(4> 

ZJN  -  Y(3)  ♦  ZJ*Y< A) 

ZEX  -  Y(7 )  ♦  ZJ*Y<8> 

ZEXP  -  Y<9)  ♦  ZJ*Y(iO> 

ZEZ  *  (ZEXP-(ZJP/VP-ZJN/VN)/EPSO/EPSO)/ZCAM 
ZHY  -  (ZJ*OMECA*EPSO*EPSC*ZEX+ZJP+ZJN)/ZGAM 
ZZ1  -  ZEZ  /  ZHY 

LOWER  GOLD  REGION 

ZKX  ■  CSGRT ( ZGAM*ZCAM+OMEGA*OMECA*ZEPSO*AMUO ) 

IF  (  AIMAG(ZKX)  LT  0.0  )  ZKX  -  -ZKX 
ZHY  ■  ( 1.  0,  0  0) 

ZEZ  -  ZKX /OMEGA/ ZEPSO 

TITANIUM.  P*  AND  P  REGIONS 

DO  40  1 DUM* 1 . 3 
I  -  7- I DUM 

IF  i  T(I)  . EG  0.  0  )  GO  TO  40 

ZKX  *  CSGRT (  ZGAM*ZGAM+OMEGA'*OMEGA*ZEPS(  I  :•  *AMUO> 


fion  n  n  n  ft  n 
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ZA  -  0.  5*  <  ZHY+Q.-iE<JA*ZEFS(  I  >*ZEZ/ZKX) 

ZB  -  0.  5*<ZHY-0MECA#ZEPS(I)*ZEZ/ZKX) 

ZHY  -  ZA#CEXP(-ZJ*ZKX*T<I> >  +  ZB*CEXP(ZJ*ZKX*T< I > > 
ZEZ  -  <ZA*CEXP(-ZJ*ZKX*T<I))  -  ZB*CEXP< ZJ*ZKX*T< I ) ) > 
1  *ZKX/0ME0A/ZEPS(I> 

40  CONTINUE 
C 

ZZ2  »  ZEZ  /  ZHY 
ZO  -  ZZ1  -  ZZ2 

C  WRITE  < 6.  200 >  ZEXI,  ZJN 

C  200 ^FORMAT  <  3X,  'ZEXI- E9.  3,  X.  E9.  3.  '  ZJN--.E9.3, 

RETURN 

C 

1000  WRITE  (6.  1010)  I PL AO 

1010  FORMAT  (  /9X.  'ERROR  (FCN2-RK.F43)  ...  IFLAC-'. 13  ) 
RETURN 
END 


IFLAC-'.  13  ) 


SUBROUTINE  ACEQN  <  X. Y. YPRIME) 

AC  EQUATIONS. 

DIMENSION  Y( 10  > . YPRIME ( 10  > 

COMPLEX  Z J.  Z 1 ,  Z2<  Z3.  Z JP  ZJN.  ZEX.  ZEXP 
COMPLEX  ZJPP.  Z JNP, ZEXPP. ZCAM. ZD 
COMMON  /BASE/PI.  DCJ.  Q. EPSO, AMUO,  OMECA 

COMMON  /8ASEM/VP.  VN, CMUN.  OMUP,  EPSC 

COMMON  /PROP /ZCAM.  ZD 

ZJ  -  <0.  0.  1  0) 

DCEX  »  Y< 1 > 

DC JN  -  Y(2> 

DCJP  -  DCJ  -  DC JN 
ZJP  “  Y<3)  ♦  Z  J*Y (4) 

ZJN  -  Y< 3 )  ♦  ZJ*Y<6) 

ZEX  -  Y<7>  +  ZJ*Y (8) 

ZEXP  •  Y<9)  ♦  ZJ*Y( 10) 

Zl  -  < Z J*0MECA+3MUN*DC JN/ VN/EPSO/EPSO ) / VN  -  GENA (DCEX) 
Z2  -  -CMUN*OC JN/VN/VP/EPSO/EPSG-CENB (DCEX ) 

Z3  -  -DC JN-CENAP ( DCEX >-DC JP*CENBP ( DCEX) 

ZJNP  -  Z1*ZJN+Z2*Z JP+Z3*ZEX+CMUN*DC JN/VN*ZEXP 

Zl  -  -<ZJ*OMECA»CMUP*DCJP/VP/EPSO/EPSG)/VP  ♦  GENB(DCEX) 
Z2  -  GMUP*DC JP/VP /VN/EPSO/EPSC+CENA ( DCEX ) 

Z3  -  DCJN*CENAP< DCEX >+DCJP*OENBP< DCEX) 

ZJPP  -  Zl*Z JP+Z2*Z JN+Z3* ZEX +GMUP*DC JP/VP* ZEXP 

Zl  -  -ZCAM*ZCAM-OMECA*OMECA*AMUO*EPSO*EPSO 
Z2  -  < Z JPP/VP-ZJNP/VN) /EPSO/EPSC 
Z3  “  ZJ*OMECA*AMUO*(ZJN+ZJP) 

ZEXPP  -  Z1#ZEX  ♦  Z2  ♦  Z3 

YPRIME < 1 )  -  <Q*CENS(X)*DCJP/VP-DCJN/VN> 

I  /EPSO/EPSC 

YPRIME (2)  -  -GENA < DCEX ) *DC JN-CENB ( DCEX > *DC JP 
YPRIMEO)  -  RE  At.  i  ZJPP  > 


•  VV  *  .*  «  vi  •  .  *  •  •  «  -  .  •  . 


*•  > 


o  o  n  non 


VPRIME<4> 
YPRIMEIS) 
YPRIME<6> 
YPR IHE<7) 
YPRIME(B) 
YPRIMEI9) 
VPRIHEUO) 

RETURN 

END 


>  AIMAG  (  Z  JPP ) 

>  REAL(ZJNP) 

>  AIMAG(ZJNP) 

'  REAL(ZEXP) 

'  AIMA3 ( ZEXP ) 

>  REALi ZEXPF ) 

-  AIMAG(ZEXPP) 


SUBROUTINE  FLDPRT  <M> 

PRINT  DC  AND  AC  FIELD  COMPONENTS  IN 
THE  ACTIVE  IMP ATT  REGION 

DIMENSION  Y< 10) , WORK ( 63 )  •  IW0RK<5) 

DIMENSION  T(6> . Z£PS(6) 

COMPLEX  ZJ. ZGAM. ZEXI. ZD. ZEPSO, ZEPS. ZKX 

COMPLEX  ZdP,  ZJN,  ZEX,  ZEXP.  ZHY,  ZEZ.  ZZ1,  ZZ2.  ZA,  ZB 

COMMON  /RKJUNK/WORK. I WORK 

COMMON  /BASE /PI.  DC  J.  0>  EPSO,  AMUO.  OMEGA 

COMMON  /BA5EM/VP. VN. OMUN. CMUP.  EPSC 

COMMON  /BOUND/XI.XF 

COMMON  /LAYER/ZEPSO, ZEPS. T 

COMMON  /PROP / ZGAM , ZD 

COMMON  /ACEXI/ZEXI 

EXTERNAL  ACEQN 

ZJ  *  (O.  0.  1.  0) 

ZEX  -  ZEXI 

UPPER  COLD  REGION 

ZKX  -  CSQRT<ZGAM»ZCAM+OHECA*OMEOA*ZEPSO*AMUO> 

IF  <  AIMAG  ( ZKX )  QT.  0.0  )  ZKX  «  -ZKX 
ZHY  «  II.  0.0  0> 

ZEZ  -  ZKX /OMEGA /ZEPSO 

TITANIUM,  N+  AND  N  REGIONS 

00  30  1-1.3 

IF  (  T<I)  .EQ  0.  0  )  GO  TO  30 

ZKX  -  CS0RT(ZCAM*ZGAM+0MEGA*0MEGA*ZEPS(I)*AMU0> 

ZA  -  0.  5* ( ZHY+OMEGA*ZEPS( I )*ZEZ/ZKX ) 

ZB  -  0.  3*(ZHY-0MEGA*ZEPS(I)*ZEZ/ZKX) 

ZHY  -  ZA*CEXP ( Z J*ZKX*T < I ) )  +  ZB*CEXP(-ZJ*ZKX*T< I ) ) 
ZEZ  -  <ZA*CEXP(ZJ*ZKX*T(I))  -  ZB*CEXP C-Z J*ZKX*T< I > > ) 
1  *ZKX/OMEGA/ZEPS< I ) 

30  CONTINUE 

ACTIVE  REGION 

NEON  -  10 
IFLAG  -  1 
RELERR  -  1.0E-8 
ABSERR  -  1.0E-12 
XI  -  XI 


uuguu  uu  uu 
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» 


ZJN  «  ZCAM*ZH Y-Z  J*OMECA*EPSO*EPSC*ZEX  ^ 

ZEXP  -  ZGAM*ZEZ-Z JN/VN/EP50/EPSC  , 

Y(l>  «  0.  0 

Y<2)  -  DC  J 

Y(3)  -00 

Y(4)  -  0.  0 

Y(S>  -  REAL (ZJN) 

Y  <  6  >  -  A I MAC (ZJN) 

Y<7)  ■  REAL ( ZEX > 

Y<8>  -  A I  MAC  <  ZEX  ) 

Y(9>  -  REAL ( ZEXP )  ' 

Y( 10)  -  AIMAC( ZEXP ) 

C 

WRITE  (1.200) 

200  FORMAT  <  //7X,  'X',9X.  'DCEX'.SX.  'DCJN'.14X.  'ACEX'. 

1  19X,  'ACJP'.  19X.  'ACJNV  > 

WRITE  (6.210) 

210  FORMAT  (  /6X,  'X',  6X.  'DCEX'.6X.  'DCJN'.llX.  'ACEX'. 

1  16X,  'ACJP '  >  1  . 

C 

WRITE  (1.220)  XI.  Y(l>.  Y(2>.  Y(7>,  Y(8>, (Y(I).  1*3.6) 

220  FORMAT  (  3X.  E9.  3.  2(2X.  E10.  4). 

1  3(2X.  E10.  4.  X>  E10.  4)  > 

WRITE  (6.230)  XI,  Y(  1 ).  Y(2>.  Y(7>,  Y(8).  Y(3>.  Y(4) 

230  FORMAT  (  X. E8. 2.  6(X,E9  3)  ) 

C 

DO  10  J-l.M  l. 

X2  -  X I ♦  J* ( XF-X I ) /M 
C 

CALL  RKF45  ( ACEQN. NEON. Y. XI . X2, RELERR. AB5ERR. 

1  I FLAG, WORK. I WORK) 

C 

WRITE  (1,220)  Xl.Y(l),  Y(2),  Y(7>,V(B>.  (Y(I).  1-3.6) 

IF  (  J.  NE.  J/10*10  .AND.  J.  NE.  M  )  CO  TO  10  “ 

WRITE  (6.230)  XI . Y( l > . Y<2> , Y<71 ,  Y(B> . Y<3> .  Y(4>  V 

C 

10  CONTINUE 
RETURN 
END 


REAL  FUNCTION  OENS(X) 

IMPURITY  DENSITY  FUNCTION  IN  THE 
IMPATT  DIODE. 


COMMON 

/ACTIVE/TN.  TP 

1 

IF  (  X 

.  OE.  0.  0  ) 

CO  TO  10 

DENS  - 

5.  0E24 

RETURN 

10 

IF  (  X 

.  CE.  TN  ) 

CO  TO  20 

DENS  - 

1.  SE23 

RETURN 

■■  ■  ■ 

L 

20 

IF  l  X 

.  CE.  TN+TP 

)  CO  TO  30 

DENS  - 

-1.  SE23 

RETURN 


L_ 


30  DENS  «  -3.  0E24 
RETURN 
END 


REAL  FUNCTION  CENA(DCEX) 

GENERATION  RATE  FOR  ELECTRONS. 

IF  (  DCEX  .  LT  JO.  >  CO  TO  10 
CENA  »  3.  3E7*£XP<-(A.  SSE7/DCEX >**2> 
RETURN 

10  GENA  ■  0.  0 
RETURN 
END 


REAL  FUNCTION  GENAP(DCEX) 

1ST  DERIVATIVE  OF  CENA(DCEX). 

IF  (  DCEX  .  LT  10.  )  CO  TO  10 

OENAP  ■  2. 0*(6.  B5E7)**2*CENA(DCEX>/DCEX#*3 

RETURN 

10  CENAP  -  0.  0 
RETURN 
END 


REAL  FUNCTION  GENB(DCEX) 

GENERATION  RATE  FOR  HOLES 

IF  (  DCEX  .LT  10.  >  00  TO  10 

OENB  -  3.  5E7*EXP<-<6.  85E7/0CEX >**2> 
RETURN 

10  GENB  -  0.  0 
RETURN 
END 


REAL  FUNCTION  GENBP (DCEX) 

1ST  DERIVATIVE  OF  GENB (DCEX). 

IF  <  DCEX  . LT  10  )  GO  TO  tO 

GENBP  -  2. 0*(6. 05E7)**2*OENB(DCEX)/DCEX**3 

RETURN 

10  GENBP  *  0.  0 
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RETURN 

END 


SUBROUTINE  ZCOMP  (USER.  Z.  EPS.  NSIC.  ITMAX.  ICOND  > 

COMPLEX  ROOT  SEARCH  :  F(Z)  «  0. 

USER  . .  .  USER  SUPPLIED  SUBROUTINE  SUBPROGRAM  IN 
THE  FORM  USER (F, Z  >  >  WHICH  RETURNS  THE 
VALUE  OF  THE  FUNCTION  F(Z>  CORRESPONDING 
TO  THE  INPUT  VALUE  OF  Z. 

Z  . .ON  INPUT,  Z  CONTAINS  AN  INITIAL  GUESS 
OF  THE  ROOT  TO  BE  FOUND. 

EPS  ...  1ST  CONVERGENCE  CRITERION. 

ABS(Z-ZP)  <  EPS  ZP  IS  PREVIOUS  Z. 
NSIC  . . .  2ND  CONVERGENCE  CRITERION. 

ABS< ( Z-ZP  > /Z )  <  10** ( — NSIG ) 

ITHAX  .  .  MAXIMUM  ALLOWABLE  NUMBER  OF  ITERATION. 
ICOND  .  .  ERROR  STATUS. 

0  -  CRITERIA  NOT  MET  (ERROR) 

1  -  1ST  CRITERION  MET 

2  -  2ND  CRITERION  MET 

COMPLEX  Z.Z1.Z2.  RT,  FO,  Fi.FP.  FRT 
IT  -  1 

DICT  -  10.  0**(-NSIG> 

FRT  -  (0.  0.0.  O) 

Z1  -  (I.  0.  1.  0) 

IF  (  CABS  ( Z  >  .EG.  0.0  )  GO  TO  10 
Z1  -  1  1*Z 
10  RT  -  Z1 
CO  TO  SO 

20  IF  <  IT  .  NE.  2  )  CO  TO  30 
FI  -  FRT 
RT  -  Z 
CO  TO  SO 
30  FO  -  FRT 

FP  -  (Fl-FO)/ ( Zl-Z ) 

Z2  -  Z  -  FO/FP 

IF  (  CABS( Z-Z2)  CE.  CABS(Zt-Z2>  )  CO  TO  40 
Zl-Z 
FI  -  FO 
40  Z  -  Z2 


CONVERGENCE  CHECK 

IF  (  CABS(Z-Zl)  LT  EPS  )  GO  TO  60 

IF  (  REAL  (  Z  )  *AIMAG  (  Z  )  .EG.  0.0)  CO  TO  43 

IF  (  ABS(REAL( Z-Z 1 ) /REAL ( Z  > )  .  LT.  DICT  AND. 

1  ABS( AIMAC ( Z— Z 1 ) /AIMAC ( Z ) )  .  LT.  DICT  > 

2  CO  TO  70 
4S  RT  -  Z 

50  CALL  USER ( FRT.  RT) 

IT  -  IT  *  1 

ITHAX  CHECK 
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